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Multiplicity distributions of intermediate mass fragments (IMF) seen in interme- 

' diate energy heavy ion collisions have been of great interest in the last ten years. The 

O ' 

distributions of single intermediate sized element are close to Poissonian and have 
been subjected to widely different physical interpretations. We calculate the multi- 
0^ \ plicity distributions in a thermodynamic model whose physical basis is equal popu- 



x 



lation in phase space which implies equilibrium statistics. We found distributions of 
single element are to a good approximation Poissonian, and such result is attributed 
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to the small emission probability of a given IMF among all the possible channels 

IT) 

' We also found significant deviation from Poisson distribution for a range of IMFs 

and suggest strong correlations as an explanation. Systematics of the multiplicity 
distributions are explored in both one component and two component models. In 
one component model, sudden changes for the multiplicity distributions are found at 
the phase transition region. In two component model, the multiplicity distributions 



have more subtleties. Results from the current model are compared with available 
data, and good agreements are found. New signals for the multi-particle correlations 
are also suggested. 

PACS numbers: 25.70.-z,25.70.Pq,25.75.Gz 

I. INTRODUCTION 

The multiplicity distributions of intermediate mass fragments (Z=3 to 20) seen in inter- 
mediate energy heavy ion collisions seem to be particularly simple. They closely resemble 
Binomial distribution. The Binomial distribution has two parameters m and p and in the 
limit m — > oo, p — > and mp — > (N) (where (N) is a constant) goes over to a Poissonian. 

n □ n 

Experimental data resemble a Poissonian |2J, |3( . This distribution has been the subject of 



widely different interpretations Q, Q, S, fl Q, , Q, Q . 

In theoretical models the multiplicity distributions can be very difficult to calculate. BUU 
(Boltzmann-Uehling-Uhlenbeck) simulations can not give accurate predictions for fragment 
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12,ll3| or Berlin model 



14, 



; lor iragmf 
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or 



24l l25| events are generated by Monte-Carlo 



productions |9j . In the Copenhagen model 1£ , 
in the lattice gas model Q H H El , 3, 22 , 
simulations and generating a multiplicity distribution would be extremely time consum- 
ing. On the other hand, new techniques have been formulated where, assuming thermody- 
namic equilibrium, various quantities of physical interest can be calculated bypassing lengthy 



Monte-Carlo simulations l2i 



3,0. 



Except for some simplifications, the physical basis of 



the model is the same as the Copenhagen or Berlin model. We therefore embarked upon 
calculating the multiplicity distributions in the thermodynamic model. 

The plan of the paper is as follows. A brief description of the canonical model is given 
in Sectjn] Sect JIIII is devoted to the multiplicity distributions of single or a range of inter- 
mediate sized clusters in one and two component model. Since decay and feeding process 
could also affect the multiplicity distributions detected in experiments, we also discussed 
them in the same section. Multiplicity distributions are strongly correlated with the system 
size, and the source size dependence is discussed in Sect JI VI The system size dependence of 
multiplicity distributions is compared with data in the same section. The multiplicity gated 
charge yield is related to the correlation (and anti-correlation) between IMFs and other 
clusters, and is discussed in Sect El The variance of multiplicity distributions is one of the 
most simple observables for cluster correlations, and is discussed in Sect 1 VII Comparison 
with data on variance is also made in the same section. In the end, we summarize the results 
in Sect JVIll The two appendices show the multiplicity tagging method for the multiplicity 
distributions, and a unified method for decay and feeding in the current model. 



II. THE THERMODYNAMIC MODEL 
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The thermodynamic model has been described in many places 
For completeness, to enumerate the parameters and the techniques needed to calculate the 
multiplicity distributions, we provide some details. We describe first the model for one 
kind of particles. Although unrealistic, the model brings out, with least complications, the 
physics of the situation readily. The generalization to two kinds is done next and follows 



the procedures given in 



29 



,3. 

If there are A identical particles of only one kind in an enclosure at temperature T, the 
partition function of the system can be written as 

Qa = ^) a , (1) 

Here u is the partition function of one particle. For a spinless particle this is uo = 
^■(27rmT) 3 / 2 ; m is the mass of the particle; V is the available volume within which each 
particle moves; A\ corrects for Gibb's paradox. If there are many species, the generalization 
is 



<k=£n^*(£ < *»'-4. ( 2 ) 



Here u>i is the partition function of a composite which has % nucleons. For a dimer % — 2, 
for a trimer i — 3 etc. Eq. (J2J) is no longer trivial to calculate. The trouble is with the 
sum in the right hand side of Eq. The sum is restrictive. We need to consider only 

those partitions of the number A which satisfy A = i n i- The number of partitions which 
satisfies the sum is enormous. We can call a given allowed partition to be a channel. The 
probability of the occurrence of a given channel P(n) = P(ni, n,2, n®....) is 

The average number of composites of i nucleons is easily seen from the above equation to 
be 

(rii) = uji^^ . (4) 



Since X^ n i = A, one readily arrives at a recursion relation 

- fc=A 

<5a = -j 2J ku k Q A _ k . (5) 
fc=i 

For one kind of particle, above is easily evaluated on a computer for A as large as 3000 
in matter of seconds. It is this recursion relation that makes the computation so easy in the 
model. Of course, once one has the partition function all relevant thermodynamic quantities 
can be computed. 
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We now need an expression for uj k which can mimic the nuclear physics situation. We 
take 

u k = ^(2vrmT) 3 / 2 x q k , (6) 

where the first part arises from the centre of mass motion of the composite which has k 
nucleons and q k is the internal partition function. For k — 1, q k — 1 and for k > 2 it is taken 
to be 

q k = exp[(W k - a(T)0 3 + T 2 k/e )/T] , (7) 

Here, as in [10"], Ho=16 MeV is the volume energy term, u{T) is a temperature dependent 
surface tension term and the last term arises from summing over excited states in the Fermi- 
gas model. The value of eo is taken to be 16 MeV. The explicit expression for cr(T) used here 
is a(T) = a [(T c 2 - T 2 )/(T 2 + T 2 )] 5 / 4 with a =18 MeV and T c = 18 MeV. In the nuclear 
case one might be tempted to interpret V of Eq. (jHJ) as simply the freeze-out volume but it 
is clearly less than that; V is the volume available to the particles for the center of mass 
motion. Assume that the only interaction between clusters is that they can not overlap one 
another then in the Van der Waals spirit we take V = Vf reeze — V ex where V ex is taken here 
to be constant and equal to Vo = A/p . 

Our objective is to calculate the probability that in an event, the probability that n 
clusters of a composite which has k nucleons is emitted. This is given by 

Pn(k)= M n Q^(^ = °) t (8) 

n! Q A 

Here QA-nkij^k = 0) is the partition function of A — nk particles but with the restriction 
that there are no composites of k nucleons. That is 

QA- nk ^k = o) = J2Il^r- ( 9 ) 

i^k % ' 

This can also be obtained by calculating the partition function for A — nk particles as before 
where all c<Vs are the same except for io k which is set to zero. 

We will also be interested in a more complicated situation. We will want to have the 
probability that n intermediate mass fragments are emitted where the intermediate mass 
fragments span a range of composites. Here we take this range to be 6 to 20. Thus n = rii 
where a is the group 6 to 20. The algebra needed to calculate this is given in Appendix 1X1 
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The decay of these clusters formed during the multifragmentation process will inevitably 
change the multiplicity distributions. Since the experiment will only observe the multiplicity 
distribution after the decay process, it is necessary to include the effects of cluster decay 
process in our considerations. The algebra is readily worked out in Appendix El The 
essential idea is to group the after-decay clusters into two groups, one with all particles that 
still stay in the group we are interested in and another with all particles that go out of 
the group we are interested in, each clusters in the two group have a modified formation 
probability. Then since we are only interested in one of the group, we will only tag the 
group of interest, and leave the other group untagged. In the end, we only need to redo the 
partition with a modified formation probability and we can find the exact solution for the 
decay process. The same idea also applies to the feeding from other clusters into the group 
of interest, with only a change of decay probability e <-» (1 — e), and a redefinition of the 
group of interest. Details of the algebra are in appendix [Bj 

These considerations can be extended to two kinds of particles, neutrons and protons 
I29I Isfl 0|. Now a composite is labelled by two indices i — > (A i: Zi), and total system is 
also labelled by two indices A — > (A, Z). The partition function for a system with Z protons 
and A — Z neutrons is given by 

® A > Z = zZIl 5 E 4 x n< - A)5(J2 Zixm-Z), (10) 

11 1 

There are two constraints: A = A^ x and Z = ^2 %i x n i • These lead to two recursion 
relations any one of which can be used. For example 

Qa,z = -g ^2 Z % UJ i Q A -A lt Z-Z l , (11) 

i 

where 

^ = ^(2™^^^$, (12) 

Here qi is the internal partition function. These could be taken from experimental binding 
energies, excited states and some model for the continuum or from the liquid drop model 
in combination with other models. The versatility of the model lies in being able to accom- 
modate any choices for qi. A choice of qi from a combination of the liquid drop model for 
binding energies and the Fermi-gas model for excited states that has been used is 

Qi = exp |- ^VoAi - a(^) 2 / 3 - k-^^ - s K - l) + — lj J . (13) 
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One readily recognizes in the parameterization above the volume term, the surface term, the 
Coulomb energy term, the symmetry energy term and contribution from the excited state. 
Coulomb repulsion between composites is sometimes taken into account by Wigner-Seitz 
approximation [lfl| . 



III. MULTIPLICITY DISTRIBUTION 



A. one component system 



The multiplicity distribution of intermediate sized clusters (6-40) is our primary interest, 
because these clusters are abundantly produced and are believed to be directly produced 



from multifragmentation process 
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Lighter clusters have 



multiple origins. Heavier clusters are often less abundant and more difficult to identify in 
experiments The essential features for multiplicity distribution is already present 

in a single component model, and we will start from this model first. We are interested in 
the multiplicity distribution of single sized clusters as well as the multiplicity distribution 
of a range of clusters. 

As is well known, the current canonical model exhibits features of a first order phase 
transition j3, 3- The divergence in the specific heat at the phase transition temperature is 
replaced by a finite peak due to finite system size effect, and the phase transition temperature 
T& is shifted to a lower value [28]. The specific heat per particle is plotted as a function 
of system size in FigCJ The temperature at which the specific heat is maximum is called 
phase transition temperature TJ, for the finite system. As will be shown later, some of the 
fragmentation characteristics change as the temperature of the fragmentation source changes 
from below phase transition temperature T < Tb to above phase transition temperature 
T > T&. Thus, the phase transition temperature provides a natural temperature scale for 
the fragmentation process in the current model. 

Fig. 12 shows the multiplicity distribution of a single cluster with mass k — 6, at three 
different temperatures T < Tb, T Tb, and T > Tb- One of the most striking features in 
Fig. 121 is that the probability of high multiplicity event is exponentially small. The Poisson 
distribution with same mean as the corresponding canonical distribution can be compared 
with the canonical results. In the linear scale used in Fig.0 it is quite difficult to see 
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FIG. 1: Specific heat Cy as a function of temperature for a system of total mass A = 200. The 
density of the fragmentation source is fixed at p = 0.05/m -3 . The peak at 6.5 MeV is defined as 
the phase transition temperature Tf, in the current one component canonical model. 

any difference between the Poisson distribution and the canonical results for multiplicity 
distribution. We have tried to fit the canonical result with either a Poisson distribution or 

r-i 

a binomial distribution as suggested by Merotto et al. |3J], the fit is always very good. The 
reason for a good fit is easy to understand. The canonical results for multiplicity distribution 
have practically only a few relevant terms, probability for all high multiplicity events are 
exponentially small. So the fitting essentially only need to fit the first few terms of the 
multiplicity distribution, and the higher multiplicity terms are just too small to affect the 
fit. In such fitting, the condition probability distribution sums up to unity is automatically 
satisfied by the fit distributions. 

In the case of T = 6.0 MeV in Fig. El we can safely neglect terms with M > 3, and 
we are left with only two significant terms in the distribution. Then it is not a surprise to 
see a good fit with Poisson distribution. Since Poisson distribution has only one parameter, 
could fit well the canonical result, we will hereafter only compare the results of the canonical 
model with the Poisson distribution. Even at T = 7 MeV, the multiplicity distribution of a 
single cluster k = 6 is well fitted by a simple Poisson distribution. And this result holds for 
effectively all small clusters as long as the system size is considerably larger than the cluster 
we are interested in. This seemingly simple distribution has been the subject of various 
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interpretations jl} 0, 0, 

In fact, the good fit with a simple Poisson distribution is not a coincidence. In a grand 
canonical model for multifragmentation process, the single cluster multiplicity distribution 
is always a Poisson distribution. The total particle number is not conserved in the grand 
canonical case, and a production of a certain cluster does not affect the yield of a second 



similar cluster at all 135 



361 ] . And each of these processes are simply combined with a factorial 



factor accounting for statistics, which then yield the Poisson distribution. For a canonical 
result to yield a Poisson distribution, we need to require small probability for single cluster 
formation, so that formation of the first cluster will not significantly affect the formation 
of a second same cluster. This condition is fulfilled when the cluster is much smaller than 
the entire fragmentation source, and there are many similar channels for particles to go into 
beside the one we are interested. In other words, the clusters in the low multiplicity events 
are effectively produced independently and little correlation exists between them. Such 
independent production mechanism together with the factorial factor for statistics always 
yields a near Poisson distribution. 
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FIG. 2: The multiplicity distribution P{M) of a single cluster with mass k = 6 as a function of 
multiplicity of this cluster M in one component canonical model. The fragmentation source has 
200 particles, and density at fragmentation is 0.05/m -3 . The phase transition temperature for this 
source is at = 6.5 MeV. The left, middle and right panels are for temperatures of T = 6, 6.5, 
and 7 MeV respectively. 
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FIG. 3: Same as El The multiplicity distribution P(M) of single clusters, k = 6, as a function of 
multiplicity M in one component canonical model, in a log plot. The left, middle and right panels 
are for temperatures of T = 6, 6.5, and 7 MeV respectively. In all three panels, the canonical results 
are shown with square symbols and the Poisson distribution with the same mean is shown with solid 
line. It is clear that canonical result for a single cluster follows closely Poisson distribution, except 
at large multiplicities. However, significant difference shows up at large multiplicities. There is a 
change of super-Poissonian distribution to sub-Poissonian distribution when temperature changes 
from below phase transition temperature to above phase transition temperature. 

Despite the good fit with a Poisson distribution, there are still significant differences in 
the canonical results and the Poisson distributions. If we look at the fine details of the 
multiplicity distribution and the Poisson distribution, we find the tails of the distributions 
differ considerably. In the log plot of the multiplicity distributions, Fig.|31 the canonical 
results show systematic deviation from the Poisson distribution. At below phase transition 
temperatures, the tail of the canonical result is significantly larger than the corresponding 
Poisson distribution with the same mean. Such a trend reverses at temperature above 
the phase transition temperature. Right at phase transition temperature, T = the 
distribution follows very closely the corresponding Poisson distribution. 

At below Tb, the system is in liquid-gas coexistence phase. The particles in the system 
could either choose to reside in the the small clusters (gas phase) or in the large clusters 
(liquid phase), and thus incur large fluctuations in the distribution. Above phase transition 
temperature, T > T&, the system is in gas phase. The particles are partitioned only in 
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small clusters. The width of the multiplicity distribution should be reduced due to finite 
system size constraint. At phase transition temperature, T = the different factors 
compete, and we are left with an almost Poisson distribution. So at the phase transition 
temperature, different constraints on the multiplicity distribution cancel, the system reaches 
a critical point. The small clusters suddenly do not feel the boundary of the system, and 
form according to the unconstrained distribution, that is, the Poisson distribution. 
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FIG. 4: The Poisson test parameter A(M) for multiplicity distribution of single clusters, k = 6, 
as a function of multiplicity M in one component canonical model. See Ea. Q14j) for the definition 
of A(M). By definition, a Poisson distribution will have A(M) = 1 for all M. The Poisson test 
parameter shows the difference in the multiplicity distribution as the temperature changes from 
below to above T^. Right at around phase transition temperature, the Poisson test parameter is 
most close to unity for a large range of multiplicities. While below and above phase transition 
temperature, the Poisson test parameter is far than from unity for most multiplicities. 

To get more quantitative description of the differences between the canonical result and 
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the Poisson distribution, it is instructive to define a Poisson test parameter 

_ (M + l)xP(M + l) 
A(M) " (M) x P(M) ' (M) 

Where P(M) is the probability of the multiplicity M event, and (M) is the mean multi- 
plicity. Such a ratio was introduced in in the discussion of multiplicity distributions. 
For Poisson distribution, this test parameter is always exactly one. For canonical results, 
this test parameter scales out the obvious power decay factor, and shows the fine details of 
deviation from Poisson distribution without reference to any artificially constructed Poisson 
distributions. 

The Poisson test parameter A(M) for the canonical results are plotted in Fig.0] With 
multiplicity distributions such as obtained for k = 6, the deviation from exact Poisson 
distribution is mostly in the high multiplicity tail. The Poisson test parameter reveals 
the quantitative deviation from Poisson distribution and such difference is readily seen at 
the high multiplicity tail. When we are dealing with a single sized cluster formed in the 
fragmentation process, P(M) reflects not only the single formation factor associated with the 
production of such cluster, it also is constrained by the weight factor for phase space of other 
clusters. With a change of M — > M + 1, as we can see from Eq.fjHJ), the first weight factor 
related to the formation of the interested cluster changes uj m /M\ — ► uj( m+1 ^ / (M + 1)!, and is 
exactly cancelled by the ratio Eq. ()14j1 ; the second weight factor related to the partition of the 
rest of the particles in the system also changes, QA-kxM^k = 0) — >• QA-kx(M+i)(&k — 0). 
When the Poisson test parameter A(M) is nearly constant for small M, we infer that there 
is a rough scaling in the second partition weight when multiplicity M changes. 

The second weight factor is just the partition of (A — k x M) particles with the condition 
of not producing the pre-selected cluster k. With the understanding that the formation of 
the pre-selected cluster is one of the many channels the system could select, ignoring such 
a channel will not significantly change the total phase space. Then it became clear that the 
scaling in QA-kxM simple suggests a system size scaling in the total partition function. At 
any given temperature, the partition pattern is not expected to change significantly when 
system size changes by a small number as long as k x M << A, so that the partition function 
will scale roughly with the system size. Such a scaling is of course going to break down when 
the size (A — k x M) does not remain large anymore, and thus finite system size effect kicks 
in. 
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While the production of a single cluster k is only a small fraction of the phase space, the 
production of a range of clusters is a much more sizable fraction of the phase space, and 
could reveal the bulk structure of the partition phase space. For single component model, 
we take the intermediate mass fragment (IMF) to be 6 < k < 40, and calculated the IMF 
multiplicity distribution, see FigsEl and El 
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FIG. 5: The multiplicity distribution P(M) of intermediate mass fragments, k = 6 — 40, as a 
function of multiplicity M in the one component canonical model. The fragmentation source has 
total mass 200, and its density is fixed at 0.05/m~ 3 . The left, middle and right panels correspond 
to temperatures below, at and above phase transition temperature. The Poisson distribution with 
the same mean is plotted together with the corresponding canonical distribution. Clearly, the 
deviation from Poisson distribution is visible even in this linear plot. The mean multiplicity for a 
range of clusters is considerably higher than that for a single cluster. 

At below phase transition temperature, T < T&, the multiplicity distribution for a range 
of IMFs is similar to the single cluster case because the the production of any IMF is small, 
and we find the Poisson fit is quite good. But when the IMFs are abundantly produced, 
T > Tb, the Poisson distribution does not describe the canonical result at all. The narrowing 
of the IMF multiplicity distribution seen at temperatures T > Tb is significant compared 
with the multiplicity distribution of a single sized cluster. Such narrow distribution reveals 
the strong correlations between the IMF clusters. Also the quick change from a fairly good 
Poisson distribution to a strong sub-Poisson distribution is a manifestation of the strong 
correlations between IMF clusters. 

The Poisson test parameter for the IMF multiplicity distribution is plotted in FigJHl Even 
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FIG. 6: The Poisson test parameter for the multiplicity distribution of a range of clusters, k=6-40, 
in the one component canonical model. The system size and breakup volume are the same as in 
FigEJ 

though the distribution is not even close to Poisson distribution at temperatures T > Tj,, 
the Poisson test parameter still reveals interesting features of the distribution. In the low 
multiplicity region, M < 10, the P(M)'s at T = Tj, is significantly lower at either T < or 
T > Tfe. Such a trend is not see in the multiplicity distribution for single IMF cluster in Fig0] 
If we look more carefully at the Poisson test parameter for M=0 in Fig0J it is clear that the 
reduction of A(M = 0) at T = 6.5 MeV is also present, but not as prominent as in Fig® 
Here, we find that a small signal in multiplicity distributions is enhanced by the selection of 
a range of IMFs, and we will also come to similar conclusions later in SectsE] and IVII For 
high multiplicities, M > 15, the lines for all three temperatures are very similar. The most 
likely clusters will be small clusters if the multiplicity is high, thus the large multiplicity 
tail of the distribution will gradually resembles the multiplicity distribution for A = 6. The 
three curves are almost parallel at very high multiplicities, which also suggests the common 
origin for the high multiplicity events. 

B. two component system 

For two component canonical model, we can also use the tagging method to find the 
multiplicity distribution for a single isotope or for a range of isotopes. The production of a 
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single light isotope is again a small fraction of the phase space, multiplicity higher than 1 
event will be only a small faction of the total fragmentation phase space. But if such rare 
event do occur, it may present different systematics due to the additional isospin degree of 
freedom. In the two component model, the partition function has two constraints, one for 
system size and another for charge, or isospin asymmetry. On the other hand, the available 
clusters has one more label for isospin too. For each fragment of given size, there are a range 
of clusters with different isospin asymmetry. 
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FIG. 7: The multiplicity distribution of Li7 and Li(all) at T = 5.2MeV, in two component model 
with A sys = 200 and Z sys = 86, in log scale. The data symbols are the results from the canonical 
model while solid lines indicate the corresponding Poisson distribution with same mean. 

The multiplicity distribution of a single isotope has similar gross features whether in 
two component system or one component system. The exponential small probability for 
high multiplicity and the almost perfect fit with a Poisson distribution are also true for 
the two component model, see Fig|7| Similar as we discussed before, such features shows 
that the production of a single isotope is just a small fraction of the total phase space. 
We can also use the Poisson test parameter to gauge the goodness of Poisson fit for the 
multiplicity distributions in two component system, and the results are shown in FigJHl The 
multiplicity distributions of single isotopes showed interesting isospin dependence, which is 
strongly related to the isospin asymmetry of the whole system. The fragmentation source 
has isospin asymmetry of 5 = (A — 2Z)/A = 0.14, the isospin asymmetry of the selected 
isotopes are <5( 5 ' 7,9 Lz) = —0.2,0.143,0.33 respectively. While the isospin of isotope 1 Li is 
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FIG. 8: The Poisson test parameter for multiplicity distribution of single clusters, Z=3, A=5,7,9, 
at three different temperatures T = 4.0, 4.6, and 5.2 MeV. The system is (A sys = 200, Z sys = 86) 
and phase transition temperature is T& = 4.6 MeV for this system. The line at A(M) = 1 is for 
exact Poisson distribution. 

most close to the whole system, the Poisson test parameter for 7 Li is much closer to unity 
than any other isotopes. The Poisson test parameters for 6 Li lie between the values for 5 Li 
and 7 Li, and is not plotted for clarity. Similarly, Poisson test parameters for 8 Li lie between 
the values for 7 Li and 9 Li. 

In the language of constrained partition problem, clusters with isospin asymmetry similar 
to the whole system, 7 Li in the current case, will feel the constraint of isospin much later 
than other clusters. For the production of 5 Li, multiple 5 'Li in the system forms an extreme 
proton-rich subsystem, which is far away from our initial neutron-rich system. To reach 
our expected isospin asymmetry in the whole system, we have to compensate this "neutron 
deficiency" with several extreme neutron-rich isotopes to get the total isospin asymmetry 
right. And we know that clusters with extreme isospin asymmetry has much lower formation 
probability as compared to other similar sized clusters. This means that 5 Li will feel the 
isospin constraint (or boundary condition) at quite low multiplicity. We may also think of 
produced isotope as building blocks from which we form the whole system. To form a system 
with isospin asymmetry 5 = 0.14, it is of course easy to use building blocks with similar 
isospin asymmetry, such as 7 Li. Then 7 Li will be more freely used than other isotopes with 
quite disparate isospin asymmetry. 

The multiplicity distributions at different temperatures also show strong isospin depen- 
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dence. As shown in FigJSl the Poisson test parameter increases with increasing temperature 
for both 5 Li and 9 Li, but the opposite is true for 7 Li. The general trend of 7 Li is in fact 
the same as the result for single component system, compare to FigHJ and should not be 
a surprise. Then why does all other isotopes will have a different temperature dependence 
for the Poisson test parameter? In the scale as used for 5 Li, we see that the lines for 7 Li 
is practically unity at all three temperatures. In this sense, we should ask why the multi- 
plicity distribution of 5 Li looks more like that of 7 Li, as temperature increases. Generally, 
isospin asymmetry in the total system is a small perturbation parameter, so that the total 
partition function may be expanded in isospin asymmetry, where the coefficient in front of 
isospin asymmetry will behave like an effective control parameter. The effect of such control 
parameter, in general, will get reduced at higher temperatures, so that the isospin difference 
between 5 Li and 7 Li became less important in the partition problem. 

In the one component model, there is a critical change from super-Poissonion to sub- 
Poissonion distributions when the temperature crosses the phase transition temperature as 
shown in Fig0] But as shown in FigJHl such a critical change is not seen in the range 
of temperatures we considered. If such a critical behavior does exist, the corresponding 
temperature should be much lower than the temperatures shown here. This is in agreement 
with conclusions in SectsfVl and IVII 



C. decay and feeding 

In the fragmentation model, clusters are produced "hot", and could decav 
particles. Such decay process will change the yield of clusters dramatically [3 
For the distribution of given clusters k, it has two effects: the cluster k could decay into 
other clusters, which we shall term as "decay" for now; and the decay of other clusters into 
cluster k, which shall be termed as "feeding". If the one of the interested cluster decays, 
we may find a multiplicity M event be identified as a multiplicity M-l event. When feeding 
happens, we will misidentify a multiplicity M event as a multiplicity M+l event. Note that 
experimentally, the incomplete coverage of the solid angles will produce the same effect as 
the decay process. 

Let us first consider the simple decay without any feeding process. For distributions 
with a mean multiplicity (M) much less than 1, the decay will keep the general shape of the 
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distribution, that is, P(M) for each higher M is order of magnitude smaller than the previous 
one. So that the difference from decay process is not readily seen in a simple multiplicity 
distribution plot. But for distributions with a higher mean multiplicity, (M) 3> 1, the decay 
process obviously changes the shape of the distribution considerably, as shown in FigEl 




FIG. 9: The effect of simple decay on the multiplicity distribution of IMFs. We have plotted the 
before and after decay multiplicity distribution of single element (Z=3), and IMFs (3 < Z < 20) 
in the left and right panel respectively. The system is (A sys = 200,Z sys = 86) at T = 5.2 MeV. 
The decay probability is assumed to be 10%. 

It is easiest to look at the changes of the mean multiplicity (M) and the variance a 2 = 
a 2 / (M) of the distribution, and the results are shown in Fig. EH The exact definitions 
for mean and variance are given in Sect JVll see Eqs|T7| and 1201 Since typical cluster decay 
probability for the produced "hot" IMF clusters is on the order of 0.10, we are assured 
that the mean value and the variance for the multiplicity distribution are not changed too 
much. But in the case of incomplete solid angle coverage, we assume the "effective" decay 
probability could be as large as 30%. 

As shown in FigJUJ the details of the multiplicity distributions are indeed significantly 
changed by the simple decay process, especially in the case of distributions of high mean 
multiplicity events. We want to see if the Poisson test parameter is changed by such decay 
process or not. For this, we have plotted the Poisson test parameter after a simple decay 
process in FigJTTl Interestingly enough, the simple decay process does not change the general 
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FIG. 10: The mean and variance of multiplicity distributions change as a function of decay prob- 
ability in two comp model. Here we have use the same multiplicity distribution of single element 
(Z=3), and IMFs (3 < Z < 20) as shown in FigUJJ but the resulting linear relation is universal. 
The mean (M) is defined in Ea. (|17l) . and scaled variance a 2 is defined in Ea.(|2U|). 

systematics of the multiplicity distribution in the Poisson test parameter plot. This suggests 
that Poisson test parameter is a good observable for comparing the current fragmentation 
model with the experiments. 

The feeding process just has the opposite effect to the decay process. In general, it will 
reduce the changes of the multiplicity distribution caused by the decay process. The decay 
probability of cluster k + 1 — ► k is generally different than the decay probability of cluster 
k — ► k — 1. For simplicity of discussion, we will assume that the decay probability of clusters 
k+1 and k are the same in the following analysis. 

The complete solution to the decay and feeding problem requires the tagging of separately 
the decay group and the feeding group, and is indeed more complicated than simple decay 
problem. Fortunately, such problem is equivalent to a new partition problem as shown in 
Appendix [Bj and could be solved exactly. Such a repartition method requires a redefinition 
of the cluster group and the cluster formation probability, but the partition part is the same 
as the tagged multiplicity distribution problem. The results for the multiplicity distribution 
after the complete decay and feeding process is shown in FigIT2l As we can see, the Poisson 
test parameter is almost not changed by the decay and feeding process. 
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FIG. 11: The Poisson test parameter A(M) for multiplicity distribution of a single element, Z=3, as 
a function of multiplicity M. The effect of simple decay will not alter the Poisson test parameter. 
The system is (A sys = 200, Z sys = 86) at T = 5.2 MeV. 

The independent production ("idp" in FiglT^j) model assumes the production of clusters 
are independent of each other. In more precise term, the probability of having N\ clusters 
in the group alpha and N 2 clusters in the group beta is: 

P(a, Aq; /3, N 2 ) « P(a, Aq)P(/3, N 2 ) . (15) 

With this combined formation probability, we can employ the same simple decay method 
as shown in the first part of Appendix [B] and do not need to build a new partition in the 
system. The independent production approximation in fact works quite well for the case of 
Z=3 cluster group, matches well to that of the exact solution up to three significant digits 
in the low multiplicity region. Actually, such statistical independence between produced 
clusters are quite expected. Since the formation of a single cluster is only a small fraction of 
the total fragmentation channel, there is little correlation between the two type of clusters, 
unless each group contains far more than one clusters, in which case, the correlation could 
build up. The deviations of the independent production model from the exact result are 
easily seen in high multiplicity events as shown in the right panel of FiglT^l 
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FIG. 12: The effects of decay and feeding on the multiplicity distribution of element Z = 3. 
The left panel shows the multiplicity distributions in linear scale, while the right panel shows the 
Poisson test parameter for the different multiplicity distributions. The line labelled "original" is 
the distribution without any consideration of either decay or feeding; the line labelled "smp decay" 
only includes the effect of decay process; the line labelled "complete" is the exact result which 
includes both decay and feeding; the line labelled "idp model" assumes the production of clusters 
are independent of each other, see Ea.(|15|). In general, inclusion of feeding process reduces the 
effect of decay and produces results more closer to "original" distribution. Such decay and feeding 
process will not change the systematics of the Poisson test parameter. Also not quite easy to see 
in the left panel, the "idp model" yields values that are up to three significant digits as good as 
the "complete" calculation in the low multiplicities M < 10. 

IV. MULTIPLICITY AND SYSTEM SIZE 



The multip icity distributions as a function of total transverse energy are measured in 
experts M- A* tb. — ene rgy is rOUg h ly COrrel ate d with the centra^ of 
the reaction, we may interpret the results as a system size dependence. But the reader is 
cautioned that when the transverse energy goes all the way down to zero, the IMF fragment 
production mechanism has changed from a central collision to a peripheral collision. To be 
consistent, we will only consider central to semi-central collision where intermediate sized 
clusters are profusely produced, and the production mechanism is better understood. The 
current thermodynamic model is, in any case, inapplicable at low temperatures where the 
fragment production is most likely sequential. 

In the one component model, the multiplicity distribution of a single cluster k = 6 is 
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shown as a function of system size in FigfHfl As the mean multiplicity for a single IMF is 
quite small, only the first few of the P(M)'s are sig nificantly above KT 2 at T = 7 MeV. 
At fixed temperature, the crossings of different P(M)'s are strongly correlated with the 
total mass of the fragmentation source. Such result is consistent with the overall features 
of multiplicity distribution for single element in We can expect each P(M)'s will rise, 
reach a maximum and then fall as system size increases. Such a trend is more clear in 
FiglT^l Here we will refrain from comparing results from one component model with the 
experimental results because the temperature correspondence are quite different in one and 
two component systems. 
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FIG. 13: Multiplicity distribution of a single cluster k = 6 as a function of system size at fixed 
temperature T = 7 MeV in one component model. 



The multiplicity distribution for a range of IMFs 6 < k < 40 is shown in FiglT4l As 
with the multiplicity distribution for single clusters, the crossing of different multiplicity 
lines for a range of clusters are also correlated with the system size. As evident in Fig lTH, 
the different multiplicity line has very similar shape for M — 10 all the way up to M = 19. 
The peak position of different multiplicity line shifts by a constant when the system size 
increase, and the peak value slightly decreases. 



Since a temperature around T 
ion reactions |3 
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— 6 MeV is often cited in intermediate energy heavy- 



44 |45j , we have plotted the multiplicity distributions 



of a single element and a range of IMFs as a function of system in FigsE3andEI]at a fixed 
temperature of T = 5.2 MeV. In producing both figures, we try to keep the system isospin 
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FIG. 14: Multiplicity distribution of a range of clusters k = 6 — 40 are plotted as a function of 
system size at fixed temperature T = 7 MeV in one component model. Only multiplicities in the 
range of 10 < M < 20 are plotted here. 
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FIG. 15: Multiplicity distribution of Z = 3 is plotted as a function of system size at fixed tempera- 
ture of 5.2 MeV in a two component model. Charge number of the fragmentation source is chosen 
so that the total system isospin asymmetry is as close to 5 = 0.1 as possible. 



asymmetry close to 5 = 0.1, but as system size and charge are integers, we can only select 
the closest integer number. As can be seen from these two figures, the variation of isospin 
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FIG. 16: Multiplicity distribution of IMFs 3 < Z < 20, is plotted as a function of system size at 
fixed temperature of 5.2 MeV in a two component model. The Charge number of the source is 
chosen so that the total system isospin asymmetry is s close to 5 = 0.1 as possible. Only lines with 
multiplicity in the range 7 < M < 17 are plotted. 

does not change much of the systematics related to system size. Other than the variation of 
isospin, we find the results for multiplicity distribution in two component system and one 
component system have similar systematics. 

At this stage, we can try to compare the systematics of multiplicity distributions with 
that from experiments. Specifically, Beaulieu et al. presented the systematics of multiplicity 
distribution of Z = 3 for the reaction Xe + Au at 50 AMeV in Fig. 4 in The transverse 
energy range plotted there is quite large < Et < lOOOMel^. But since the low transverse 
energy data are most probably produced in peripheral reactions where the cluster production 
mechanism may be of a sequential nature, we will only concentrate on the more central 
collisions where the transverse energy is in the range of 500MeV < Et < lOOOMeV. The 
total system size could be assumed to be proportional to the total transverse energy, while 
the system temperature is roughly constant, thus we estimate the system size to be 120 < 
A sys < 240. This assumption is of course too simplified. There may also be a change of 
temperature when the transverse energy Et changes. For simplicity, we also assumed the 
total system isospin asymmetry is roughly 5 = 0.10. 

As can be seen from Fig. 4 in Q , the mean multiplicity of Li is considerably lower than that 
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in FiglTBI To reach lower mean multiplicity, we either have to choose a smaller system or a 
lower temperature. On top of these two uncertainties, we do have another freedom to change 
the decay probability for Li element to simulate the incomplete coverage of emitting angles. 
So we choose to fix the system size as we have estimated, and try to fit the experimental 
result with only variation of temperature and decay probability. Additionally, a large decay 
probability and a lower temperature have similar effect in reducing the mean multiplicity, 
with the temperature effect being the more prominent of the two. So we may constrain 
ourselves to a high temperature and lower decay probability. Since in fragmentation models, 
a temperature of around 5 MeV is often used 0, l^l] , we will try to find a sensible decay 
rate while keeping temperature as close to 5 MeV as possible. With the current published 
data, we have found no reference to get an estimate on the effective decay probability, we 
will just assume it is between 15% and 30%. Since this decay probability also includes the 
incomplete coverage of the detector system, this range is not too far away from reality. 
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FIG. 17: Multiplicity distribution, for single element Z = 3, is plotted as a function of system size 
in a two component model. We have assumed 15% decay probability for all Li isotopes and the 
temperature of the fragmentation source is fixed at T=4.6 MeV. Such result roughly corresponds 
to the multiplicity distribution in the range 500MeV < Et < lOOOMeV at the left panel of fig. 4 
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We have plotted the multiplicity distribution of Li as a function of system size in FigfTTl 
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FIG. 18: The Poisson test parameter, for single element Z = 3, is plotted as a function of system 
size. We have assumed 15% decay probability for all Li isotopes and the temperature of the 
fragmentation source is fixed at T=4.6 MeV. This result roughly corresponds to the right panel 
of fig.4 in with the range 500MeV < Et < WOOMeV. Note, our definition of Poisson test 
parameter A(M) scaled out the obvious factor of mean multiplicity, see Ea.(|14j). 

A temperature of T = 4.6 MeV and a decay probability of 15% seem to give a good fit to the 
data in Fig.4 in Combinations of either a slightly lower temperature and lower decay 
probability or slightly higher temperature and higher decay probability also gives similar 
results. 

The Poisson test parameter for the same distribution in FigH3is plotted in FiglTHl The 
this result could be compared with the right panel of Fig.4 in Note in the definition of 
Poisson test parameter Eq. (jl4|) . we have scaled out the mean multiplicity. The results seem 
to agree very well with the systematics in 

The multiplicity distribution of a range of IMFs, 3 < Z < 20, is also presented in |3j. In 
the reaction of Ar + Au at 110 AMeV, multiplicity distributions of IMFs are measured in 
the transverse energy range < Et < HOOMeV^. As before, we will try to reproduce the 
multiplicity distributions in the more central collisions with 550 < Et < HOOMeV", and the 
system size is estimated as 110 < A sys < 220. Again, we have assumed the total system 
isospin asymmetry to be around 5 = 0.10. 

As shown in FigEH if we assume a temperature of 3.5 MeV and a decay probability of 
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25%, the systematics of the multiplicity distribution reasonably produces the result in jjj. 
Here the temperature of the source is significantly lower than the previous fitted 4.6 MeV 
for Xe + Au reaction. 
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FIG. 19: Multiplicity distribution, for a range of IMFs 3 < Z < 20, is plotted as a function of 
system size for two component system. We have assumed 25% decay probability for all Li isotopes 
and the temperature of the fragmentation source is fixed at T=3.5 MeV. This result roughly 
corresponds to the range of 550MeV < Et < HOOMeV in the top panel of Fig. 3 in 



In FigEDl we have plotted the Poisson test parameter for the same multiplicity distri- 
butions in Fig|19l This result can be compared with that in FigJTHl In the single element 
case, the Poisson parameters are always quite close to unity; while its magnitude changes 
considerably for a range of IMFs. This is also a manifestation of the correlations between 
IMF clusters, which will give rise to stronger signals for a large range of IMFs than for a 
single IMF (see also discussion in Sect. IV1J) . The line for A(M = 0) shows oscillations as a 
function of system size. This oscillations may be related to the "long range correlations" in 
the cluster production, where particular patterns of cluster formation matches well with the 
total system size and thus get enhanced at the optimum system size, and may get suppressed 
when mismatch happens. We will see more of this matching effect in Sect. 
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FIG. 20: The Poisson test parameter, for a range of IMFs 3 < Z < 20, is plotted as a function of 
system size for two component system. The IMF selection and Et correspondence are the same 
as in Fig|19l We have assumed 25% decay probability for all Li isotopes and the temperature of 
the fragmentation source is fixed at T=3.5 MeV. 

V. MULTIPLICITY GATED YIELD RATIO 



The most elementary measurements are inclusive cross sections. These are given by the 
yield (n k ), see Eq.(dJ). The isotope yield ratios are often used in heavy-ion reactions because 
some experimental uncertainties are cancelled in the ratios. Examples include single iso- 
tope yield ratios, double isotope vie 



reactions \3l 



34, 



38, 



41 



43, 



44, 45, 



djratios, and ratios of same isotope but from different 
It is also possible to make more complicated mea- 
surements: measure (n^) subject to the condition that the event has a given IMF multiplicity 
M. We can define then a ratio: 



R k {M) 



(n k (M)) 
(n k ) 



(16) 



where (n k (M)) is the yield when the IMF multiplicity is M and (n k ) is the usual (ungated) 
yield. 

In general, if we fix the IMF multiplicity to be lower than the mean IMF multiplicity, then 
the IMF production is artificially suppressed. If any other clusters are positively correlated 



28 




FIG. 21: The multiplicity gated to ungated particle yield ratio Rk(M), for a single component 
system at around the phase transition temperature region. The three panels from left to right are 
for the same system of A sys = 200 particles at different temperatures T = 5, 6 and 7 MeV, while 
the phase transition temperature TJ, = 6.5MeV. For IMF in the range of 6 < k < 40, the mean 
multiplicities at the three temperatures are 0.12, 1.55, and 12.72 respectively. 

with the IMF clusters, the yield of such clusters will also get suppressed. On the other 
hand, if the selected multiplicity is higher than the mean multiplicity, the IMF clusters get 
enhanced, and so will the positively correlated clusters. 

The IMF multiplicity gated and ungated yield ratio Rk(M) for a system of 200 particles 
at three different temperatures are plotted in Fig|21l The three panels from left to right are 
results at temperatures well below phase transition temperature T = 5MeV, close to phase 
transition temperature T = 6 MeV, and above phase transition temperature T = 7MeV. 
The corresponding mean multiplicity for the systems from left to right are (M) = 0.12, 1.55, 
and 12.72 respectively. We can see a change of shape for different multiplicity lines when 
temperature changes from below to above phase transition temperature. Such a shape 
change is an exhibition of the critical behavior in the single component system. 

Aside from the obvious shape change in Fig. |^ we can also find the correlation between 
IMF clusters and heavier mass clusters. At 5 MeV, most of the heavier mass clusters are 
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FIG. 22: The multiplicity gated to ungated element yield ratio Rk(M), for a single component 
system, with A sys = 200 and T = 7MeV. 

positively correlated with the IMFs when M > 0, and the line for M = is strongly 
suppressed; while at 7 MeV, each of the P(M)'s has an positively correlated region and a 
strongly reduced correlation at heavier masses. The only exception is the multiplicity M = 
line, where only positive correlations is observed. 

As we have indicated before, the selection of a range of IMF clusters will typically enhance 
the signals. Here we may test the enhancement by changing the selection for IMF clusters 
and see if there is any change in the pattern of yield ratio. We varied the IMF range and 
plotted the lines for yield ratio in Figl2*2l The multiplicity M = lines seem to be specially 
interesting, aside from simple change of magnitude in the yield ratio, certain oscillating 
structure shows up when the range of IMF is large enough. The appearance of oscillating 
structure is due to the "long range multiparticle correlations" between heavier clusters and 
will be elucidated in more detail here. 

In the case of high temperature T = 7 MeV> the system is predominately in the 
gas phase, and small clusters are profusely produced, and these small clusters are highly 
correlated. All other larger clusters will have negative correlation with the IMFs, and a 
select of M = for IMFs suppress the productions of IMFs and enhance the production of 
larger clusters. The additional oscillations in the cluster yield ratio, represent the multiply 
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correlations between clusters of quite different size, and thus are of "long range" nature. For 
example, the line labelled "6 — 60" in Fig|22l has a significant dip at around k = 115 while 
the line "6-40" does not show a dip at the same fragment size. This dip must come from the 
exclusion of IMFs in the range 40 < k < 60. This could happen when the source prefers to 
fragment into two large pieces, one of size 115 and another of size between 40 and 60, and 
also some small clusters. Since the probability of M = event is rare, it is not a surprise to 
see the exotic fragmentation process. 
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FIG. 23: The multiplicity gated to ungated charge yield ratio Rz(M) in a two component model. 
The total system is A sys = 200, Z sys = 86 at T=4 MeV. IMFs are defined as clusters with 
3 < Z < 20, and the mean IMF multiplicity is (M) = 6.45 in this system. 

The gated to ungated charge yield ratio in a two component model shows similar features 
as in the one component model, but the change of shape does not show up at temperatures 
close to the phase transition temperature. Figl2"3l shows the typical yield ratio for a two com- 
ponent model, where the system size is taken as A sys = 200, Z sys = 86, and the temperature 
of the fragmentation source is well below phase transition temperature T = 4MeV < T&. 
And if the shape change exists in the two component model, the corresponding temperature 
should be well below temperatures we are interested here. This result is consistent with the 
conclusion of Sect J VII 

Another interesting feature is the appearance of multiple oscillation in the line M = 
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FIG. 24: Same at in Fig. I2.3I The multiplicity gated to ungated charge yield ratio in a two 
component model where the IMF range selection is varied. The system is A sys = 200, Z sys = 86 
at T=4 MeV. The different lines correspond to different IMF range higher cut-off from Z < 10 to 
Z < 40, while the lower cut-off is fixed at Z > 3. 

in FigJ2Hl The appearance of two dips at Z — 43, 66 are pointing to two different exotic 
fragmentation channels, one is related to two clusters of size Z ~ 20 and Z ~ 66, and 
another is related to three clusters of size Z ~ 43 and two of Z ~ 20. The significant 
charge difference between such clusters could be labelled as "long range", because most 
of the correlations we considered are between similar sized clusters, or "short ranged". It 
should be noted that the probability for these exotic fragmentation channels is quite small 
as compared with that of regular channels which also produce the heavy cluster Z = 66 
together with some IMF clusters. If we allow the production of small clusters in the range 
of 6 < Z < 10, the dip at Z = 66 for will not show up, because the most of the channels 
producing a Z = 66 are accompanied by several small clusters. On the other hand, the line 
for M = 1 also shows oscillations and the phase of oscillations is just opposite to the line of 
M = 0. The difference in the phase of the two lines is another manifestation of the strong 
correlation in the exotic fragmentation channels. 

Fig]2*4l shows the effect of different IMF range selections on the yield ratios for fixed IMF 
multiplicity M = 0. As with one component model, the enhanced yield of heavy elements 
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are due to the negative correlations between IMFs and heavy clusters. The oscillation of 
the M = lines are more prominent in the two component model. The strong oscillation 
and the strong dip in the M = line suggest strong correlations between large cluster and 
small cluster. If we look at the largest charge value for the dip to happen Z max (dip), we find 
this values for different lines in FigE] are of equal distance. The position of the dip is at 
Z max (dip) = 74, 66, 56, and 46 respectively, and adding the cut-off value for the IMF range 
selection, we find the sum to be 84, 86, 86, 86. These exotic fragmentation channels are very 
similar to asymmetric binary fission. 



VI. VARIANCE AND PHASE TRANSITION 



The variance is one of the most simple cumulants for the distributions, and is often used 
as a synonym for fluctuation or correlation in the system. Higher order cumulants were also 



studied before 



47, 



49, 



50|. Coinciden tally , the variances for charge fluctuations 



tanv 
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transverse momentum fluctuations (see |5_3j and reference therein) are also proposed for 
relativistic heavy-ion collisions in the study of a possible phase transition to the quark-gluon 
plasma phase. Experimentally, the variance is easiest to measure. The variance is related to 
the correlations between particles, and the appearance of a maximum variance is sometimes 
interpreted as a critical process in intermediate energy heavy-ion collisions 

Q, Q- 

The strong correlations between clusters, or fluctuations, are sometimes attributed to the 
phase transition in the fragmentation source 47. |. In this section, we will test the relation 
between correlations and phase transition and try to fit the fragment correlation results from 
experiments. 

The mean multiplicity and the variance are defined as 

<M> = ^MP(M), 

M 

a 2 = ^M 2 P(M) - (M) 2 , 



(17) 



M 



where M is the multiplicity of the interested particles. Note that from the definition of P(M), 
the definitions above are complete equivalent to the usual definitions for mean multiplicity 
and variance: 



(M) = £» 
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o 2 = (5>) 2 >-<I» 2 , (18) 

For Poisson distribution, the variance is exactly equal to the mean value, so it is instructive 
to use the scaled variance o 2 = a 2 / (M), for which a Poisson distribution just gives unity. 
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FIG. 25: Specific heat Cy and scaled variance a 2 as a function of temperature in one component 
canonical model. The temperature at maximum specific heat is defined as the phase transition 
temperature in the canonical model. For a system of size A sys = 200 and the density p = 
0.05/m -3 , the phase transition temperature is at = 6.5 MeV. The peak for maximum IMF 
correlations is at T = 6.4 MeV. IMFs are defined as clusters of size 6 < k < 40. 



In the one component model, occurrence of maximum scaled variance (or simply max- 
imum correlation) is simply related to the phase transition. As seen in Fig OBl the phase 
transition temperature and the maximum variance almost occur at the same temperature. 
We can also see the close relation between temperatures for maximum correlation T c and for 
maximum specific heat T), in FigQol The difference between the two temperatures becomes 
smaller when the system size increases, and we could expect the two temperatures be the 
same for a infinite large system. Similar result for the peak and valley structure for variance 
has been reported for a lattice gas model [54 1 . 

As already observed in Sect llll Al the multiplicity distributions shows systematic changes 
when the temperature of the system changes from below to above phase transition temper- 
ature. Plotted in FigE3 is the scaled variance for multiplicity distributions of IMFs. The 
line labelled with "6-6" only selects clusters of size k = 6, and the corresponding variance 
does not show a strong variation across the critical temperature region. Such a result is just 
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FIG. 26: The temperatures for max specific heat and max correlation are plotted as a function of 
system size in the one component model. The density of the source is at p/po = 1/3. 

as expected from our previous analysis on the multiplicity distributions of single clusters, 
where a simple Poisson distribution fits well the multiplicity distribution. The change of 
the variance for a single cluster is quite small, and poses difficulties when comparing with 
experiments. However, as we select larger range of clusters, the variance shows significant 
structure of maximum and minimum in the region of phase transition temperature. The 
line with IMF selection of 6 < k < 40 shows a strong peak at close to the phase transition 
temperature, and a quick drop right above phase transition temperature. Such a behavior 
is an interesting feature, and is much easier to compare with experiments. The chan ge o f 
variance across the phase transition region has been observed before in a similar model pfj] . 

The buildup of the variance is due to the correlations between clusters. If the clusters 
are uncorrelated, we have: 

o 2 = <(£^) 2 >-<Q>>)> 2 , 

i i 
* id 

= £°? = I>K> ( 19 ) 

i i 

where the cross terms are uncorrelated and cancel out. The scaled variance is 

rmn(7f) < 7 = = < ™ax(7f) , (20) 
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Thus, the sum of uncorrelated distributions could not yield variance higher than the largest 
variance of the individual distributions, nor could it yield variance lower than the smallest 
variance of the individual distributions. Any significant enhancement or reduction of the 
variance as observed in FigEZJ must come from correlations between clusters. 
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FIG. 27: The scaled variance is plotted as a function of temperature in the phase transition region. 
Different lines correspond to different IMF selections. The change of the variance for a single 
cluster k = 6 is quite small across critical temperature, while the change is significantly enhanced 
when including a large range of clusters. 



As can be seen from Figs 123 and |2H1 the variance peak is a continuous function of system 
size and temperature, so that it forms a ridge in the contour plot in Fig|2Hl The left 
hand side of the correlation ridge has variance always above one, and the right hand side 
are always below one, with a significant minimum at above but still close to the ridge. If 
experimentally we find the variance are below one and changes continuously with system size, 
the only allowable temperature in the one component model is above the phase transition 
temperature. If the variance value is larger than 0.2, which is the value for the minimum in 
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FIG. 28: The scaled variance as a function of temperature and system size for the one component 
canonical model. The shade of darkness represents the value of the scaled variance and the corre- 
spondence is shown by the insert scale. The contours regions are drawn in steps of 0.05 in scaled 
variance, while contour lines are drawn in steps of 0.15 starting from a 2 = 0.275. The contour 
region corresponding to a 2 = 0.6 is darkened while that for a 2 = 0.45 is hatched. This result 
may be compared with Fig. 3 in Q] in the transverse energy range 700MeV < Er < lAOOMeV. 

this model, then there is practically only a very narrow range of allowable temperature just 
above the phase transition temperature. 

We may compare the variance result with that from an experiment, Fig.3 in [lj. The 
region we are interested in are the more central collisions, with 700MeV < Et < 1400MeV , 
in Xe + Au reactions at 110 AMeV. The corresponding system size is estimated to be 
150 < A sys < 300 assuming the linear relation between transverse energy and system size. 
The scaled variance changes from a 2 = 0.6 to a 2 = 0.45 in the data. The contour regions 
corresponding to these value of scaled variance are indicated in FigEHl By projection of the 
hatched region to the temperature axis, we find only a narrow range of allowable temperature 
for the system. Such a tight constraint on the temperature is a unique feature of the one 
component model. 
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FIG. 29: Contour plot of the scaled variance as a function of temperature and charge of the 
source Z sys , where system size is fixed A sys = 110. The dark colored band indicates the region 
with a 2 = 0.82, corresponding to the variance measured by Moretto et al, in Ar+Au at 110 
MeV/nuleon. The hashed band indicates the region with a 2 = 0.75, which is the decay corrected 
value for the variance. Each contour band corresponds to a change in variance by 0.035, while 
contour lines are also drawn in steps of 0.105 starting at a 2 = 0.35. The temperature obtained 
here is around 3 MeV, if the isospin asymmetry of the source is around 5 = 0.15. This temperature 
is slightly lower than 3.5 MeV assumed from the multiplicity distribution result in Fig ll9l in Sect llVI 

For two component model, we have one additional degree of freedom in the isospin di- 
rection and situation became more complicated. The parameter space we are interested in 
is {A, T, 5} here. We need to understand the variance dependence on isospin asymmetry 
5, or in other word, system charge Z sys . Figsl2T?l and 1301 show the variance of IMFs for a 
system of A sys = 110 and 220 particles respectively. As can be seen, the isospin dependence 
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FIG. 30: Contour plot of the scaled variance as a function of temperature T and charge of the 

source Z sys , where system size is fixed A sys = 220. The dark colored band indicates the region 

~ Q 
with a 2 = 0.62, corresponding to the variance measured by in Fig. 4 in |3], in Ar + Au at 110 

MeV/nuleon. The hashed band indicates the region with a 2 = 0.5, which is the decay corrected 

value for the variance. Each contour band corresponds to a change in variance by 0.02, while 

contour lines are drawn in steps of 0.06 starting from a 2 = 0.31. The temperature we find here is 

around 3.5 MeV if the isospin asymmetry of the source is estimated as 5 = 0.15. This temperature 

is in agreement with that from multiplicity distribution results in Sect llVI 

of the variance is quite smaller in the higher temperature region. And there is a significant 
variance peak at low temperature of around T = 2MeV. It is generally believed that at very 
low temperature, a sequential emission scenario is more appropriate than the simultaneous 



fragmentation model 



12| . So we will try to look for regions where the temperature is not 



too low, generally T > 3 MeV. 
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We may compare the variance results with data from experiment. We will concentrate 
on the variance data as show in Fig.l in |2J. Again, we will only consider more central 
collisions and assume simple linear relation between transverse energy and system size, 
specifically in the region of 550 < Et < HOOMeV^, or system size of about 110 to 220. The 
variance changes from a 2 = 0.82 to a 2 = 0.62, and the corresponding contour region are 
indicated with dark gray color in Figs I2H1 and EH In the discussion in Sect. IIV[ we have 
assumed the a decay probability of 25% to account for simple decay of IMFs and incomplete 
coverage in the detection of IMFs. As shown in FigJTUl the variance changes linearly with 
the decay probability, so we can make simple corrections for the decay process. The decay 
corrected variances are a 2 =0.75 and a 2 = 0.5. The contour regions corresponding 
to the decay corrected variances are shown with hatched patterns in FigsEHl and EH In 
the smaller system, a correction of decay only slightly increases the corresponding value 
for temperature, but the increase is more significant in the larger system. The extracted 
temperatures from FigsEHl and EDI depend on the isospin asymmetry of the system, and 
we find around T = 3MeV and T = 3.5MeV for the two system if isospin asymmetry 
of the system is around 5 = 0.15. In the smaller system, a significant change in isospin 
asymmetry of the system would yield considerably lower temperature, which we have tried 
to avoid. In the larger system, A sys = 220, the temperature we may get from pure decay 
corrected variance is in the range of 3 < T < 3.8 MeV if the isospin of the system is in 
the range of 0.07 < 5 < 0.20. The temperatures extracted for the larger system seem to 
agree well with the fitted temperature T = 3.5MeV in Fig El But the smaller system 
does show some deviation. This discrepancy maybe due to the simple assumption of linear 
relation between transverse energy and system size, or the unrealistic assumption of constant 
temperature when transverse energy changes, or a combination of the two. As we have 
mentioned in Sect II VI an increase system size, an increase in temperature and a reduction of 
decay probability have similar effects. In light of the isospin effect on variance in Figll9| we 
might expect a slightly larger system size and a temperature of around 3 MeV to better fit 
both the multiplicity distribution result and the variance result at around Et = 550 MeV. 
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VII. SUMMARY 

In a canonical multifragmentation model, we used a recursive procedure to solve exactly 
the constrained partition problem. A recursive method for the multiplicity tagged partition 
function is derived. Multiplicity distributions of single or a range of clusters are calculated 
exactly. 

The Multiplicity distributions of single intermediate mass fragment (IMF) are found to 
be well approximated by a Poisson distribution. Details of the multiplicity distribution of 
single clusters show a critical behavior near the phase transition temperature in the one 
component canonical model. Following the suggestion of we defined a Poisson test 
parameter for the multiplicity distributions. This Poisson test parameter reveals details of 
the fragment production process near the phase transition. In a two component model, the 
multiplicity distribution of single clusters shows additional sensitivity to isospin asymmetry 
of the fragmentation source. The selection of a large range of clusters introduces strong 
correlations between clusters, thus the critical behavior shows up more strongly in the one 
component model. When comparing with experiments, the decay and feeding effects will 
significantly change the multiplicity distributions of single clusters or a range of clusters. 
However, the change of the mean multiplicity and the variance changes linearly with the 
decay probability. The production of single clusters is mostly uncorrelated, so that the 
multiplicity distribution of two clusters is well approximated by the simple product of each 
multiplicity distributions. 

The multiplicity distributions of single and a range of clusters have a smooth dependence 
on the system size. The crossings of different multiplicities have a strong correlation with 
system size if we assume the temperature is fixed. The multiplicity distribution of a range 
of clusters seems to show repeated patterns as a function of system size, which is very 
suggestive of repetition of similar fragmentation patterns as the system size increases. We 
can fit the data from Q, [| with simple assumption about the system size in a two component 
canonical model. An increase of system size, increase of temperature, and decrease of decay 
probability have similar effect on the multiplicity distribution. By requiring the system size 
fixed, we can find a single temperature that fits the data well. The Poisson test parameter 
for the distributions shows smooth dependence on system size. The Poisson test parameter 
for single clusters is always quite close to one, while the selection of a large range of clusters 
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will generate more appreciable differences in the Poisson test parameter. 

The IMF multiplicity gated charge yield shows a change shape in the one component 
model, while such shape change did not show up in the two component model. The IMF 
multiplicity gated charge yield also shows correlations between IMF and other heavy clusters. 
Oscillating structure appeared in the zero IMF multiplicity gated yields, which is related to 
the rare process of binary fission type of fragmentation process, where the two fragments 
are strongly correlated. 

The variance is one of the most simple cumulants of the multiplicity distributions, and 
is also a simple measure of the correlations between IMFs. In a one component model, the 
critical behavior of the variance appears in the phase transition region. The selection of a 
large range of IMFs is shown to significantly enhance the change of variance across the phase 
transition region. The two component model is used to explain the correlation result from 
experiment |2|. The fitted parameters of temperature and decay probability are consistent 
with that from fitting the multiplicity distributions. 

Overall, we believe the multiplicity distributions and related signals provide new insights 
into the multifragmentation process in heavy-ion reactions. 
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APPENDIX A: MULTIPLICITY TAGGING 

Let us denote each species that can emerge from the breakup as a,b,c,d,... etc. Let us 
divide these into two groups: a and (3. A given species belongs to either a or (3 but not to 
both. It is obvious that the partition function of the whole system which has A nucleons 
can be written as 

Qa= Y,Qa 1 {oc)Q A2 {(5)5{A 1 + A 2 -A). (Al) 

Ai,A 2 

Here Qa x (o) means the partition function of A\ nucleons which can be partitioned into 
clusters in the a group. For example, a could comprise all the intermediate mass fragments. 



42 



We want to do a multiplicity tagging on the partition within the a group. For this, we 
define 



Qa 1 (a, M) = II 5 E in * - A ^ 5 E n < - M ) » (A2) 



iga i 

then clearly 



Q A (a) =^Q A (a,M). (A3) 



,1/ 

Starting from Eq. (jA2|) and taking advantage of the identity 



M 

tea 

we can build up a recursive formula for Qa{&, M). 

Q A (a,M) = Ef^EM n^ 5 E^-^ 5 (E^- M ) 

iGa \ fc€a / i 

''^ - L / - L rij! 

i, tea i^tfe 



m E ^rdbw n 



X 



(Tlk — 1J! , n i ! 
i, tea j^fc 

<5( J] mi + A;(n fc - 1) - (A - k)) S(J2 *k + K - 1) - (M - 1)) 

= ^E Wfc ^- fc ( a ' M_1 )- ( A5 ) 

tea 

Going back to the full system (a + (3), the relevant probability is given by 

P(a,M) = -J- QA 1 (a,M)Q A2 ((3)5(A 1 + A 2 -A). (A6) 

APPENDIX B: DECAY AND FEEDING 

We may start from the simple case of decay only problem. To make it even simpler, 
we will assume the clusters in the group a have equal probability e to decay into clusters 
outside the group a. The decay of clusters within the group a is not considered since they 
do not change the multiplicity distribution. 

If before decay, the multiplicity distribution is P(a, N), then after independent decay of 
each of these cluster in the group a, the multiplicity N event will become multiplicity M 
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event (M = 0, 1, 2, N), with a weight of 

P(N — > M) — CZe N ~ M (l - e) M , (Bl) 
Then the probability of multiplicity M event after decay is 

Pd(a, N — »■ M) = ^ P(a, N)P(N — ► M) , (B2) 

TV 

However, when the decay probability of the clusters in group a is not uniform, the above 
equation will not work. For this more general case, we have to tag the multiplicity of each 
of the clusters in the group a, which will be a formidable task for even a moderate sized 
system. 

Because of decay, if the number of composite of type i in group a was this can turn 
into Ni composites in group a with probability 

P{n . ^ N .) = C^er^l - Ci)^ • (B3) 

The total after decay tagged partition function is: 

Q >,M) = £n^(".^)n^ 

i j i 

The total weight factor for cluster % can be simplified: 



71A Hi 



6iUJ, 



\rii-Ni 



[(1 - 



(B5) 



Here we may view the first term as a partition weight for a new type of cluster i', where 
the formation probability is (b\ = (e^), N- — m — N^, and the second term in the partition 
could be viewed as a weight for the after decay cluster i, with formation probability uj-i = 
[(1 - edfa)]. With this, 

We will denote the after decay clusters % E a as the stay group a s , and the after decay 
clusters i' as a new group a D . 
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With the help of Eq. (jB6|) . the total tagged partition Eq. (jB4|) is simplified into 



X 



Nil 11 M! f 1 i\T,-! 
5(Y, Ni - M) <f(£ iNi + £ i'iVj + £ jiV, - A) (B7) 



Where the new formation probability is defined as 



UJ; 



\ 



(1 - 6i)ui, i G a s ; 

e^, i e a D ; (B8) 



Since we are only interested in the group a s , where tagging is needed, the rest of the clusters 
could be combined into a new group called (3' = a D + other. 

l£a a i£a a i£a° 

x n -jJt 6 ^ jNj - M) 5{Ai +a *- a ^ 

je? r jep' 

= Qm (<* S , M)Q A2 {{3') 5(A 1 +A 2 — A). (B9) 

Ai,A 2 

And the multiplicity distribution is now 

P d^ m) = 1_ Q Al (a s ,M)Q A2 (P')6(A 1 + A 2 ~A) (BIO) 

WA A l> A2 

Notice Eq. ()B10|) resembles Eq. (jA6|) . the differences are in the grouping of clusters, and 
in the formation probability (u instead of uj). 

The decay of clusters in group (3 could also feed into the group a, and therefore changes 
the multiplicity distribution of after-decay multiplicity distribution of the clusters in group 
a. Since each of these clusters decays independently with a different decay probability, we 
have to tag each of these before-decay clusters in group (3 as well as a. The multiplicity for 
the cluster (j,rij) will be multiplied by the decay weight factor 

P( nj - = C%ef{l - e,) n ^ . (Bll) 

Notice Eq. (|B11J) differs from Eq. ()B4j) by only the exchange of e <-» (1 — e). 
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The total after-decay tagged partition function is 



<#-(<,, m) = y. n - *<> n ^ ^ n 



j'e/3 J k<£a,k£f3 

x5(j2 w + E j N i + J2 kNk ~ A ) ^E ^ + E N i - M ) ' ( B12 ) 

i j k i j 

As can be readily seen, the factorization of Eq. (|B6)1 works for both i 6 a and j G (3. We 
only need a redefinition of the formation probability o> 



C^i = < 



1 - ei)u)i, i G a b ; 

eiUi, i G a D ; 

[1 - ei)ui, % G (3 s ; 

em, i G (3 D ; 

u u i <£a,i <£ (3. 



(B13) 



And the tagged after-decay partition function is 



Q r(a, M) = sn^n ^(E + e - *o 



ke/3 1 



TT fa)"' 



n 



X 



TT ^^TT 

J-J- JV.-! 11 iVj! iV m ! 

5(J2 + E ^ + E fciVfc + E m + E mN ™ - A ) ( B14 ) 

i j k I m 

Define the new set of interest a' = a s + (3 D , and irrelevant (3' = a D + (3 s + other. Then we 
have 
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